Comparison of approaches for predicting solute transport: sandbox experiments     

Illman WA  Berg SJ  Yeh TC 《Ground water》2012,50(3):421-431

The main purpose of this paper was to compare three approaches for predicting solute transport. The approaches include: (1) an effective parameter/macrodispersion approach (Gelhar and Axness 1983); (2) a heterogeneous approach using ordinary kriging based on core samples; and (3) a heterogeneous approach based on hydraulic tomography. We conducted our comparison in a heterogeneous sandbox aquifer. The aquifer was first characterized by taking 48 core samples to obtain local-scale hydraulic conductivity (K). The spatial statistics of these K values were then used to calculate the effective parameters. These K values and their statistics were also used for kriging to obtain a heterogeneous K field. In parallel, we performed a hydraulic tomography survey using hydraulic tests conducted in a dipole fashion with the drawdown data analyzed using the sequential successive linear estimator code (Yeh and Liu 2000) to obtain a K distribution (or K tomogram). The effective parameters and the heterogeneous K fields from kriging and hydraulic tomography were used in forward simulations of a dipole conservative tracer test. The simulated and observed breakthrough curves and their temporal moments were compared. Results show an improvement in predictions of drawdown behavior and tracer transport when the K tomogram from hydraulic tomography was used. This suggests that the high-resolution prediction of solute transport is possible without collecting a large number of small-scale samples to estimate flow and transport properties that are costly to obtain at the field scale.  相似文献   

An analytical solution for heterogeneous and anisotropic anticline reservoirs under well injection     

Hund-Der Yeh  Chia-Chen Kuo 《Advances in water resources》2010

This study develops a mathematical model for describing the steady-state head response to fluid injection into a fully penetrating well in a heterogeneous and anisotropic anticline reservoir. In the model, the upper boundary of the anticline reservoir is approximated by a form of step change in reservoir thickness and the domain of the reservoir is divided into two regions with different hydraulic conductivities. By virtue of the properties of Fourier series, the method of separation of variables is employed to develop the analytical solution of the model.  相似文献   

A new analytical solution solved by triple series equations method for constant-head tests in confined aquifers     

Ya-Chi Chang  Hund-Der Yeh 《Advances in water resources》2010

The constant-head pumping tests are usually employed to determine the aquifer parameters and they can be performed in fully or partially penetrating wells. Generally, the Dirichlet condition is prescribed along the well screen and the Neumann type no-flow condition is specified over the unscreened part of the test well. The mathematical model describing the aquifer response to a constant-head test performed in a fully penetrating well can be easily solved by the conventional integral transform technique under the uniform Dirichlet-type condition along the rim of wellbore. However, the boundary condition for a test well with partial penetration should be considered as a mixed-type condition. This mixed boundary value problem in a confined aquifer system of infinite radial extent and finite vertical extent is solved by the Laplace and finite Fourier transforms in conjunction with the triple series equations method. This approach provides analytical results for the drawdown in a partially penetrating well for arbitrary location of the well screen in a finite thickness aquifer. The semi-analytical solutions are particularly useful for the practical applications from the computational point of view.  相似文献   

Two-dimensional stochastic analysis of flow in leaky confined aquifers subject to spatial and periodic leakage     

Hund-Der Yeh  Ching-Min Chang   《Advances in water resources》2009,32(11):1601-1608

This paper addresses the issue of flow in heterogeneous leaky confined aquifers subject to leakage. The leakage into the confined aquifer is driven by spatial and periodic fluctuations of water table in an overlying phreatic aquifer. The introduction of leakage leads to non-uniformity in the mean head gradient and results in nonstationarity in hydraulic head and velocity fields. Therefore, a nonstationary spectral approach based on Fourier–Stieltjes representations for the perturbed quantities is adopted to account for the spatial variability of nonstationary head fields. Closed-form expressions for the variances of hydraulic head and specific discharge are developed in terms of statistical properties of hydraulic parameters. The results indicate that the spatiotemporal variations in leakage leads to enhanced variability of the hydraulic head and of the specific discharge, which increase with distance from any arbitrary reference point. The coefficient of leakage and the spatial structure of log transmissivity field and of the amplitude of water table fluctuation are critical in quantifying the variability of the hydraulic head and of the specific discharge.  相似文献   

Basin-scale transmissivity and storativity estimation using hydraulic tomography     

Kuhlman KL  Hinnell AC  Mishra PK  Yeh TC 《Ground water》2008,46(5):706-715

While tomographic inversion has been successfully applied to laboratory- and field-scale tests, here we address the new issue of scale that arises when extending the method to a basin. Specifically, we apply the hydraulic tomography (HT) concept to jointly interpret four multiwell aquifer tests in a synthetic basin to illustrate the superiority of this approach to a more traditional Theis analysis of the same tests. Transmissivity and storativity are estimated for each element of a regional numerical model using the geostatistically based sequential successive linear estimator (SSLE) inverse solution method. We find that HT inversion is an effective strategy for incorporating data from potentially disparate aquifer tests into a basin-wide aquifer property estimate. The robustness of the SSLE algorithm is investigated by considering the effects of noisy observations, changing the variance of the true aquifer parameters, and supplying incorrect initial and boundary conditions to the inverse model. Ground water flow velocities and total confined storage are used as metrics to compare true and estimated parameter fields; they quantify the effectiveness of HT and SSLE compared to a Theis solution methodology. We discuss alternative software that can be used for implementing tomography inversion.  相似文献   

Integration of principal components analysis and cellular automata for spatial decisionmaking and urban simulation   总被引：2，自引：0，他引：2  

黎夏  叶嘉安 《中国科学D辑(英文版)》2002,45(6):521-529

This paper discusses the issues about the correlation of spatial variables during spatial decisionmaking using multicriteria evaluation (MCE) and cellular automata (CA). The correlation of spatial variables can cause the malfunction of MCE. In urban simulation, spatial factors often exhibit a high degree of correlation which is considered as an undesirable property for MCE. This study uses principal components analysis (PCA) to remove data redundancy among a large set of spatial variables and determine ‘ideal points’ for land development. PCA is integrated with cellular automata and geographical information systems (GIS) for the simulation of idealized urban forms for planning purposes.  相似文献   

Uncertainty in applying the linear poroelasticity model to field situations as a result of periodic loading in heterogeneous aquifers          下载免费PDF全文

Ching‐Min Chang  Hund‐Der Yeh 《水文研究》2015,29(10):2454-2462

The solid Earth's surface frequently experience changes in total stresses as a result of periodic loading. When the fluid‐saturated porous media deform in response to changes in stress, the induced variations in pore volume affect the pore water pressure. The fluid flow therefore occurs in response to the gradient in the induced excess pore water pressure. This work aims at quantifying the spatial variability in excess pressure head produced by the periodic loading accounting for the variation of log hydraulic conductivity (lnK). It is important for the rational management of groundwater resources. A closed‐form expression is developed by the nonstationary spectral approach to analyse the influence of the statistical properties of lnK process, the hydraulic parameters, and the spatial position. The general stochastic framework outlined in this work provides a basis for assessing the impact of statistical properties of input aquifer parameters on the output variability (or uncertainty). Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Embedding isochronous cells overland flow module into MODFLOW     

Bo Zhang  Mei Hong  Zuhao Zou  Tian‐Chyi Jim Yeh  Pin Cheng 《水文研究》2013,27(26):3833-3841

This study integrates a simple overland flow module (isochronous cells model) with the river module of MODFLOW such that temporal and spatial interactions between stream flow and groundwater can be simulated using net rainfall data of a watershed. The isochronous cells model is an efficient travel time runoff approach based on geographic information system (GIS) that considers both spatial and temporal variations of net rainfall through hill slope of the watershed. This overland module is easily coupled with MODFLOW river routing module. Specifically, the stream flow from the isochronous cells model is directly assigned to both sides of river cells of the MODFLOW model. Such an integration of MODFLOW and isochronous cells model is especially useful in watersheds where river flow data are limited. The feasibility of this integrated model was demonstrated using a case study in the middle and downstream regions of the Yitong River watershed, China. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Soil water balance model for precipitation-induced shallow landslides   总被引：3，自引：2，他引：1  

Hsin-Fu Yeh  Cheng-Haw Lee 《Environmental Earth Sciences》2013,70(6):2691-2701

Precipitation infiltration is one of the most significant triggering factors for slope failure occurrence in many places around the world. Knowledge of the mechanisms leading to precipitation-induced slope failures is of great importance to the management of landslide hazard. In this study, a soil water balance model is developed to estimate soil water flux during the process of infiltration from rainfall data, with consideration of storm periods and non-storm periods. Two important assumptions in this study are given: (1) instantaneous uniform distribution of the degree of saturation and (2) a linear relationship between evapotranspiration and the related degree of saturation. For storm periods, the Brooks and Corey model estimates both the soil water retention curve and soil water parameters. The infiltration partition is employed by an infinite-series solution of Philip in conjunction with the time compression approximation. For none-storm periods, evapotranspiration can be derived for the moisture depletion of soil water. This study presents a procedure for calculating the safety factor for an unsaturated slope suffering from precipitation infiltration. The process of infiltration into a slope due to rainfall and its effect on soil slope behavior are examined using modified Mohr–Coulomb failure criterion in conjunction with a soil water balance model. The results indicate that the matric suction, which is closely related to slope stability, is affected by the degree of saturation controlled by rainfall events.  相似文献   

红树林湿地植被生物量的雷达遥感估算   总被引：19，自引：0，他引：19  

黎夏  叶嘉安  王树功  刘凯  刘小平  钱峻屏  陈晓越  何执兼  覃朝锋 《遥感学报》2006,10(3):387-396

根据雷达后向散射系数建立了红树林湿地植被生物量的估算模型,并运用遗传算法确定其中非线性模型的最优参数.对比分析表明,雷达后向散射系数模型比NDVI模型在植被生物量估算中有更高的精度.使用NDVI指数有可能导致某些植被类型的生物量估算出现较大的误差.这是因为一些具有密集冠层的草本植被（例如互花米草等）有比红树林高得多的NDVI值.而雷达遥感所具有的侧视特点及一定的穿透能力能有效地获取植被的垂直信息,大大减低植被生物量估算的误差.  相似文献